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ABSTRACT 

We investigate the possibility of describing the spiral pattern of the 
Milky Way in terms of a model of superposition 2- and 4-armed wave 
harmonics (the simplest description, besides pure modes). Two com- 
plementary methods are used: a study of stellar kinematics, and direct 
tracing of positions of spiral arms. In the first method, the parameters of 
the galactic rotation curve and the free parameters of the spiral density 
waves were obtained from Cepheid kinematics, under different assump- 
tions. To turn visible the structure corresponding to these models, we 
computed the evolution of an ensemble of A^-particles, simulating the 
ISM clouds, in the perturbed galactic gravitational field. In the sec- 
ond method, we present a new analysis of the longitude-velocity (/ — v) 
diagram of the sample of galactic HII regions, converting positions of 
spiral arms in the galactic plane into locii of these arms in the / — v dia- 
gram. Both methods indicate that the "self-sustained" model, in which 
the 2- armed and 4-armed mode have different pitch angles (6° and 12°, 
respectively) is a good description of the disk structure. An important 
conclusion is that the Sun happens to be practically at the corotation 



circle. As an additional result of our study, we propose an independent 
test for localization of the corotation circle in a spiral galaxy: a gap 
in the radial distribution of interstellar gas has to be observed in the 
corotation region. 

Subject headings: Galaxy: kinematics and dynamics - Galaxy: structure 
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1. Introduction 

A good understanding of the large-scale 
spiral structure of the Galaxy has not been 
reached up to the present. Georgelin & Georgelin 
(1976), herafter GG, derived a 4- armed pat- 
tern, based on an analysis of the distribution 
of giant H 11 regions. According to Vallee 
(1995) most researches support the 4-armed 
pattern, although there are discordant opin- 
ions; for instance Bash (1981) finds the pat- 
tern to be 2-armed, similar to the first spi- 
ral structure model proposed by Lin and Shu 
(1964). But even if we accept that there are 
strong observational evidences in favor of a 4- 
arms structure, theoretical and observational 
difficulties remain, and the observation of ex- 
ternal galaxies strongly recommends that we 
look for more complex solutions, which should 
include arms with different pitch angles. 

It is now accepted that the distance to the 
center is about 7.5 kpc (eg. Oiling & Mer- 
refield, 1998), which implies a smaller rota- 
tion velocity of the Local Standard of Rest 
than previously estimated, and as a conse- 
quence, the rotation curve decreases beyond 
the solar radius (Amaral et al., 1996, Honma 
& Kan-Ya, 1998). From the rotation curve, 
the epicyclic frequency and the radius of the 
inner and outer Lindblad resonances can be 
derived. Amaral and Lepine (1997, here- 
after AL) concluded that the 2-armed pattern 
could exist between 2.5 and 12 kpc, which is 
about the range of the observed spiral pat- 
tern, but a 4-armed pattern could only exist 
between about 6 to 11 kpc. AL proposed a 
representation of the galactic spiral structure 
superposition of 2- and 4-armed wave 
patterns, to account for the number of arms 
and for the existence of arms over a wider 
range of radius. In AL's model, two arms of 
the 4-armed component are coincident with 



the 2-armed component, so that the Galaxy 
looks 4-armed. However, there arc theoret- 
ical arguments discussed in the present pa- 
per showing that the pitch angle for the 4- 
armed pattern should be different from that 
of the 2-armed component. A Fourier analysis 
of external galaxies in terms of spiral modes, 
performed by Puerari & Dottori (1992) in- 
dicates as well that in the cases in which the 
2-armed and 4-armed patterns are prominent, 
they indeed have different pitch angles. These 
arguments suggest that at least in a range of 
radius, the Galaxy might present a structure 
like 2-arms plus 4-arms with different pitch 
angles, which could look like 6 arms. An at- 
tractive aspect of such a model with arms of 
different pitch angles is that it naturally ac- 
counts for bifurcations of arms. In external 
galaxies branching of arms is a widespread 
phenomenon, e.g. the galaxy M 101, where 
there are both branching arms and bridges 
between them. The structure of the Milky 
Way is perhaps similarly complicated. An- 
other fact that stimulates us to look for more 
refined solutions to the spiral structure is that 
there are observed arms in the Galaxy that 
do not fit in a simple 4-arms structure. For 
instance, outside the solar circle, around lon- 
gitudes 210° - 260°, three arms are clearly vis- 
ible, in HI and in IRAS sources (Kerr, 1969, 
Wouterloot et al., 1990), but they are not 
expected from a 4-arms model that fits the 
tangential directions of the inner parts of the 
Galaxy (e.g. Ortiz & Lepine, 1993). 

The main goal of this paper is to inves- 
tigate if a superposition of 2-|-4-armed wave 
harmonics is a good representation of the spi- 
ral structure of the Galaxy. Two different ap- 
proaches are used. In one approach, we anal- 
yse the Cepheid kinematics, using a model 
that takes into consideration the perturba- 



3 



tion of the stellar velocity field by the spi- 
ral arms, to derive the structural parameters 
for the gravitational field of the galactic disk. 
And, since the structure of the perturbation 
potential does not necessarily correspond to 
the visible structure, we performed a many- 
particle simulation of gas dynamics in this po- 
tential, to predict the visible structure. In 
the second approach, we directly analyse the 
visible structure by means of the sample of 
HII regions, since these objects are recognized 
as the best large-scale tracers of the galactic 
structure. We present a new analysis of the 
observed longitude- velocity (/ — v) diagram of 
the best up-to-date sample of H II regions, de- 
riving from it the position of the spiral arms, 
and we compare the observed / — v diagram 
with the theoretical ones, computed from our 
gas dynamics simulations. 

The Cepheid kinematics analysis is per- 
formed for two different models of 2-arms -|- 4- 
arms structures, and the results are compared 
with previous calculations of pure 2-arms and 
pure 4-arms. An important result is that for 
the more reallistic model, the Sun lies very 
close to the corotation region. Furthermore, 
our simulation revealed the interesting effect 
of gas pumping out from the corotation. This 
suggests that in spiral galaxies, there must 
be a gas deficiency in a region near the coro- 
tation circle. This phenomenon explains the 
lack of atomic hydrogen in a ring-like region 
near the solar circle derived by Kerr (1969) 
and by Burton (1976), which remained unun- 
derstood up to now. We propose this effect 
as an independent test for localization of a 
corotation circle in external spiral galaxies. 



2. Method of estimation of structural 
parameters 

The derivation of the spiral wave param- 
eters is based on the statistical analysis of 
stellar motion in the Galaxy (sec Crczc & 
Mcnnessier, 1973, Mishurov & Zcnina, 1999, 
hereafter MZ, Mishurov et al., 1997, hereafter 
MZDMR and papers cited therein). Note 
that we look for the parameters of the struc- 
ture of the galactic gravitational field, which 
is not directly visible. However, the gravita- 
tional field determines the stellar motion, and 
in particular, the spiral perturbations of the 
field deviate the stellar motion from rotation 
symmetry. Hence, analysing the stellar veloc- 
ity field in the framework of a perturbation 
model, we can derive both the parameters of 
the density waves and those of the galactic 
rotation. 

Unlike the above cited papers, let us rep- 
resent the gravitational potential (po of the 
Galaxy as the sum: 

<PG^<Po + <PS2 + VS4 (1) 

where (po is, as usual, the unperturbed ax- 
isymmetric part of the potential that de- 
termines the Galaxy equilibrium whole 
{d(po/dR = fl^R, fl{R) is the angular rotation 
velocity of galactic disk, R is the galactocen- 
tric distance) and (fsm {m — 2 or m — A) are 
the m—th harmonics of the perturbations due 
to spiral density waves. According to Lin et 
al (1969): 

(pSm^ -^mCOsiXm), (2) 

where Am is the amplitude of the m — th har- 
monic, Xm is the wave phase: 

Xm^m [cot{im) ln{R/RQ) -'&]+ Xem , (3) 

m is the azimuthal wave number, i.e. the 
number af arms for a given harmonic, is 
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the corresponding pitch angle of the arms, R, 
are the polar coordinates with the origin at 
the Galactic center, Rq is the solar galacto- 
centric distance, Xem is the initial phase or 
the wave phase at the Sun position. Prom 
Eq.(3) it is seen that this last value fixes the 
m — th harmonic position relative to the Sun. 

In accordance with Eq (1) we can write 
the perturbed stellar velocity for the radial 
(directed along galactocentric radius) vr and 
the azimuthal components as follows: 

Vr = /i?2 cos(x2) + fR4 cos(x4) , (4) 

Vi) = U2 sin(x2) + /w sin(x4) , (5) 

where Jr^ and f^rn are the amplitudes of 
m — th harmonic. These quantities are re- 
lated to the parameters of the density waves 
by some formulas (see Lin et al, 1969). 

Substituting Eqs (4-5) into Eqs.(l,2) of MZ 
and using the statistical method described 
by MZ and MZDMR, we derive the parame- 
ters of the rotation curve (fi©. A, Rq^q ) and 
the parameters of the spiral waves {fRm, f^mi 
im, X&m) over the observed stellar velocity 
field. Then by means of the density wave the- 
ory (Lin et al, 1969) we compute the differ- 
ence Afl„i between the angular rotation ve- 
locity oi m — th pattern and the rota- 
tion velocity of the Galaxy at the Sun po- 
sition {AQjn = ^pm ~ ^o)- By equating 
^{Rcm) = ^pm we find the displacement ARm 
of the Sun relative to the corotation radius 
Rem {ARm = Rem — Rq, for Computational 
details see MZ and MZDMR). 

Notice here some features of our task. It 
is widely believed that the spiral structure of 
our Galaxy is generated by a bar in the galac- 
tic center (e.g., Marochnik & Suchkov, 1984). 
So, the angular rotation velocity of the pat- 
tern is determined by the bar rotation. Hence, 



the value does not depend on m, and the 
corotation radius does not depend on m ei- 
ther. 

Further, the nature of the perturbed grav- 
itational field could be two-fold. First, in Lin 
et al.(1969) theory, the spiral waves are self- 
sustained and the disturbances of the gravi- 
rational field are mainly due to the den- 
sity waves which propagate in the galactic 
disk. The corresponding dispersion relation 
imposes a connection between 12 and i^. In- 
deed, the solution of the dispersion relation 
for spiral waves relative to the radial wave 
number k.^ in the vicinity of corotation does 
not depend on m (Shu, 1970, Mark, 1976). 
Since cot (im) = kmR/iTi, we have 

cot(i2) — 2cot(i4) . (6) 

Therefore, the two-armed pattern is tighter 
wound (about twice for small pitch angles) 
than the four-armed one. 

The above argument is strictly held in the 
vicinity of the corotation circle (according to 
Creze & Mennessier, 1973; MZ and MZDMR; 
AL, etc. the Sun is situated just near the 
corotation). However, since the radial wave 
number and consequently the pitch angle are 
slowly varying functions of R (Lin et al, 1969) 
the relation (6) can be used for a sufficiently 
wide region around the corotation radius. 

In the second approach the perturbations 
of the gravitational potential are mainly due 
to a bar in the galactic center. According to 
the models in which 2 of the 4-arms compo- 
nents coincide with the 2-arms components 
(eg. Englmaier & Gerhard 1999, Al): 

X4 = 2x2; H = i2- (7) 

In what follows we shall analyze both these 
approaches. We call them the self-sustained 
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model (approach 1) and the bar-dominated 
model (approach 2). 

Another important peculiarity of our task 
is that some of the wave parameters: /^2! f&2, 
fm-ifm-iii {h is fixed automatically in the 
both models), Xq2 and XqAi obey nonlinear 
statistics. In the self-sustained model Xqi ^^'^ 
Xo4 are considered to be independent quan- 
tities; in the bar-dominated approach, these 
values are connected by Eq (7). But if we fix 
the pitch angles and the initial wave phases, 
the task becomes a linear one over other the 
parameters. So, the strategy to localize the 
global minimum for the residual (5^) in the 
second approach is just the same presented 
by MZ and MZDMR. However, in the self- 
sustained model it shghtly changes. In this 
case we have to fix 3 independent quantities: 
^2, X&2 and Xoi) and look for the minimum of 
residual over other the parameters by means 
of the least square method (we denote this 
minimum by A = mini^^^^^^^^^ ((52)). Then 
we change values of ^2, X&2 and Xqa and again 
derive A, and so on. After that, we construct 
the net A as a function of ^2, X02 and X04. Of 
course, we cannot imagine this function vi- 
sually, but we can construct the surface, say, 
A(x02, X04) for a set of values 12 and check by 
eye localization of the global minimum of the 
residual. After the minimum is localized, by 
means of a linearization procedure (Draper 
& Smith, 1981) we define more exactly the 
searched parameters and the covariation ma- 
trix of errors (see details in MZ). 

3. Results of Cepheid kinematics anal- 
ysis 

The above described procedure was ap- 
plied to a sample of Cepheids which represent 
the best observational material for solving the 
above formulated problem. For the Cepheids 



both the line-of-sight velocities (Pont et al, 
1994, Gorynya et al, 1996, Caldwell k Coul- 
son, 1987) and proper motions (HIPPAR- 
COS, ESA, 1997) are available. Let us discuss 
the results derived for the above approaches 
separately. 

Self- sustained model. For ilhistration in 
Fig.l the surfaces ^{xq2,Xq4) are given for 
3 values of i2-0ne can easily see the global 
minimum in the vicinity of 12 ~ —6°. The fi- 
nal values of the parameters with their errors 
are given in Table 1. 

First of all, we notice the significant de- 
crease of the residual min{5'^) in this case in 
comparison with the one for the single har- 
monic m — 2 or m — A (see Table 1 of MZ, 
runs no 3 and 8). As it was shown in the 
cited paper, the inclusion of the spiral per- 
turbation in stellar motion happens to be sig- 
nificant. However the authors could not make 
a choice between the two alternatives: pure 2- 
or 4-armed patterns. 

Now by means of F-test we can show that 
the previous hypotheses that the pattern can 
be represented by only one harmonic (m = 2 
or m = 4), should be rejected in favor of the 
alternative hypothesis that the pattern is well 
represented by superposition of 2-|-4-armed 
pattern. In other words, the representation 
of the galactic structure by a superposition of 
2-|-4-harmonics is clearly preferable to the one 
that uses a single harmonic. 

The quantities Afl^n and Ai?^ can both 
be calculated from the parameters obtained 
for m = 2 and for m = 4. It is not hith- 
erto obvious that for different m the quan- 
tities occur to be the same as we have sup- 
posed above. But our calculations lead to 
very close values for the corresponding quan- 
tities: Afl2 = 0.15 km s~^kpc~^ and Af24 = 
Q.l'ikm s~^kpc~^] accordingly Ai?2==-0.03 /cpc 
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and Ai?4=-0.04 kpc, the standard (i.e. 68%) 
confidence intervals being for AQ2 '■ —0.61 to 
1.02 km s-^kpc-^; for Ail^ : 0.13 to 0.2Akm 
s-^kpc-^; for AR2 : -0.21 to 0.13 kpc; for 
Ai?4 : —0.05 to 0.03 kpc (the confidence in- 
tervals were estimated by means of numerical 
experiments described by MZ). Hence, in the 
model under consideration the Sun is practi- 
cally situated at the corotation circle, slightly 
beyond it. 

We cannot estimate with any reasonable 
accuracy the values of the amplitudes of the 
spiral gravitational field A2 and A4 (this is 
expected from linear perturbation theory). 
However, their ratio is derived very precisely: 
A2/A4 = 0.79, the standard confidence inter- 
val being 0.77 to 0.80. 

The locus of minima for Lpsm arc shown in 
Fig. 2; they are the lines of constant phase 
Xm on the galactic plane corresponding to 
minipsm- From this figure the pattern may 
be thought to be 6-armed one. However, this 
is not the case in our model, in which the po- 
tential perturbation is represented by a cosine 
functions. Indeed, simple computation shows 
that for the above derived parameters the sum 
VS2 + ¥^54 has at most 4 minima over i? for 
a fixed R. It would only be possible to ob- 
tain 6 minima (in the frame of a 2+4 armed 
model) if the potential were represented by 
some function presenting sharp minima, con- 
trary to the cosine function. The visible struc- 
ture derived by means of particle-cloud sim- 
ulations and given in Sec. 4 supports this 
point of view. Of course, the actual structure 
of the Galaxy may occur to be more compli- 
cated e.g., due to higher wave harmonics or to 
special effects at the corotation (Mark, 1976). 
However, these possibilities are beyond the 
present investigation. 

Bar- dominated model. In this case the 



results are quite different. The pattern ro- 
tation velocities are: for m = 2 Qp2 — 
35.0 kms^^kpc^^ and for m = 4 Qp4 = 
29.2 kms~^kpc~^ . So that the main require- 
ment of the model (independence of pattern 
rotation velocity and the corotation radius 
from m) is not held. Indeed, for parame- 
ters of rotation curve of Tables 1 and the 
above value for ^2^2 the corotation radius i?c2 
does not exist as a real number at all (this is 
mainly because the second derivative Rq^Iq 
is negative). On the other hand for m = 4 
i?c4 = 5.4 kpc. Further, in this case the ratio 
A2/A4 ^ 8.21 is very large. So, the visible 
pattern happens to be 2-armed! 

At last, the residual 5"^ ^ 220 is signifi- 
cantly greater than in the previous approach 
and occurs to be very close to the values for 
pure m = 2 and pure m = 4 solutions (see 
Table 1 of MZ 1, runs no 3 and 8). Hence, it 
is impossible to make a choice between pure 
171 — 2, pure m = 4 or superposition of 2-1-4- 
modes in the bar-dominated approach. 

The above statistical analysis of the large- 
scale stellar kinematics of Cepheid stars. leads 
us to the conclusion that the preferable solu- 
tion for the spiral structure of the Galaxy is a 
superposition of self-sustained 2+4:-\mTmonics 
of density waves, and that the Sun is situated 
very close to the corotation circle. 

4. Visible leirge-scale structure of the 
GalELxy 

In the previous Section we derived the 
structural parameters for the gravitational 
field of the Galaxy. However, as it was men- 
tioned above, it is not directly seen. To make 
visible the above derived structure, the evo- 
lution of a gas cloud ensemble in the galac- 
tic gravitational field pertubed by spiral arms 
will be next considered. Following Roberts & 
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Hausman (1984), we simulate the interstellar 
clouds by ballistic particles moving in a given 
gravitational field with a potential ipc defined 
by Eqs (1-3). Let us briefly describe the for- 
mulation of the task. 

At the initial moment of time {t — 0) the 
spiral perturbations are assumed to be absent 
{^Sm = 0). N particles {N = 2-10^) are 
uniformly distributed over a disk within R < 
13 kpc. Each particle is given the local ro- 
tation velocity, disturbed by a chaotic veloc- 
ity with one-dimensional dispersion 8 km s~^. 
For t > the spiral perturbation is " switched 
on". Our task is to compute the reaction of 
the system on this perturbation (the A^-body 
problem for particles moving in an external 
field). 

The parameters both for unperturbed and 
for perturbed potential were taken from Table 
1. Since the derived rotation curve is valid 
in a restricted region, for i? > 9.4 kpc we 
continue it by a fiat part. 

The self-sustained galactic waves are well 
known to exist between the inner and outer 
Lindblad resonances. Since we do not take 
into consideration a bar in the galactic cen- 
ter, the spiral perturbation was cut off for 
R < 2 kpc. For the spiral gravitational am- 
plitude A2 we assumed the "standard" value: 
2A2Cot{t2)/n\^R^Q = 0.05 (Lin et al 1969). 
In the simulation, collisions between the in- 
terstellar clouds occur, the collisions being 
energy-dissipating. For the computation of 
this process, we used the method described in 
detail by Roberts & Hausman (1984). For the 
cloud cross-section we adopted a typical value 
for the H I clouds. Mutual clouds gravitation, 
or other effects like interaction of the clouds 
with expanding envelopes of supernovae, etc., 
were not taken into account. We also restrict 
ourselves to consideration of two-dimensional 



particle motion in the galactic plane. All com- 
putations are performed in a frame of refer- 
ence corotating with spiral arms. 

The result of our simulation of particle- 
cloud dynamics in the spiral gravitational 
field for the superposion of 2-1-4 self-sustained 
density wave harmonics is shown in Fig. 3. 
There is good agreement in major features 
with Fig. 3 of Efremov (1998). In a signif- 
icant part of galactocentric distances the pat- 
tern looks hke a 4-armed one. But we do not 
face with the problem of too short arms, as it 
would be in the case of pure m = 4 harmonic 
(see AL). Our pattern reflects the complicated 
picture often observed in external galaxies, 
e.g. arm bifurcation or their overlaping. 

5. Gap in the galactic gaseous disk as 
an indicator for the corotation cir- 
cle 

One of the most important conclusion of 
Sec. 3 is that the Sun lies very close to the 
corotation radius. In this Section we present 
a new test which turns possible to localize di- 
rectly the position of the corotation circle in 
a spiral galaxy. 

Many years ago Kerr (1969) payed atten- 
sion to a ring-like region which is markedly 
deficient in neutral hydrogen, with radius 
slightly greater than the solar distance from 
the galactic center (see also Simonson, 1970). 
This result was later supported in more de- 
tail by Burton (1976). He showed that there 
is a very clear gap in radial distribution of 
atomic hydrogen in our Galaxy at i? ~ 11 kpc, 
whereas in the old scale, used in that paper, 
i?0 = 10 kpc (see Fig. 6 of Burton, 1976). In 
general, the gap reminds the Cassini gap in 
Saturnian ring. 

It is natural to connect this gap in the ra- 
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dial ISM distribution with the process occur- 
ing at corotation: the gas is pumped out from 
the corotation under the influence of the grav- 
itational field of spiral arms (see also Suchkov, 
1978). The quahtative explanations is as fol- 
lows. It is well known, that when the gas 
flows through the galactic density waves, a 
shock arises in the medium (Roberts, 1969; 
Roberts & Hausman, 1984). Since the galac- 
tic disk rotates differentially and for R < Rc 
Q > Qp, the gas in this region overtakes the 
spiral wave, entering it from the inner side. 
In the shock the clouds are decelerated, and 
fall towards the galactic center. For R > Rc 
Q < Qp, and the process is inverse. Here 
the wave overtakes the gas and pushes it. 
So, the clouds pass to an orbit more remote 
from the galactic center (see also Goldreich 
& Tremaine, 1978 and Gor'kavyi & Fridman, 
1994). 

Our simulation of gas cloud dynamics in 
spiral gravitational field directly demonstrates 
this phenomenon. We show in Fig.4 the radial 
gas distribution (< n > is the particle concen- 
tration averaged over a circle) for t = and 
t — 3.0 (the time is given in rotation period at 
the solar distance) for the best parameters of 
Sec 3. The gap in the ISM distribution at the 
corotation radius is well seen. Comparison of 
Fig.4 with Fig. 6 of Burton (1976) shows a 
close similarity between them. 

This result enables us to explain another 
problem as well. It is well known that the ro- 
tation velocity of the disk presents a sharp 
minimum near the solar galactocentric dis- 
tance; the minimum appears independently 
of the tracer being gaseous (eg. Honma & 
Kan-ya, 1998) or stellar (Amaral et al, 1996). 
This phenomenon, could be thought to be 
understood in terms of the velocity pertur- 
bation from the galactic density waves, or 



from the rising effect of a dark matter com- 
ponent at distances larger than the solar ra- 
dius. Amaral et al. exclude these hypotheses, 
in their discussion of the nature of this min- 
imum. Now, if there is a ring-like region de- 
void of gas, in principle the rotation velocity 
could not be measured in that region, using a 
gas tracer. Similarly, if one selects short-lived 
stars as tracers, these stars are not expected 
to form and to exist inside that region, be- 
cause of the gas deficiency and of the very 
low velocity of the gas with respect to the 
spiral pattern, which turns the gas compres- 
sion (the star-formation process) inefficient. 
It is therefore probable that the gas clouds 
and the stars that we observe inside the gap 
are objects with non-circular orbits that in- 
vade the gap; they are observed close to their 
maximum elongation and therefore, present 
smaller velocity than the circular one, in the 
direction of rotation. 

So, a ring-like gap in the galactic gaseous 
disk, and possibly also a sharp minimum in 
the rotation curve, may serve as an indepen- 
dent indicator for localization of the corota- 
tion circle in a spiral galaxy. 

6. New spiral structure of the Galeixy 
derived from H II data 

The H II regions are the best tracers of 
the large-scale spiral structure, since they can 
be observed at large distances, and unlike H 
I, they are sharply concentrated in the arms 
(e.g. GG). We performed a new analysis of 
the sample of H II regions of the Galaxy, with- 
out making use of the results of the theoretical 
models discussed in this paper, except for one 
hypothesis, namely, that the structure can be 
represented by a superposition of 2-|-4-armed 
patterns. Therefore, the results of this section 
constitute an independent test of the previous 
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ones. 

The procedure adopted is to trace spiral 
arms in the galactic plane, and to transform 
X positions along the arms into locus of 
the arms in the / — v diagram, by means of 
the rotation curve. By varying the parame- 
ters of the arms, we looked for the best fit to 
the / — V diagram. It is well known that the 
arms that are situated inside the solar circle, 
transform into narrow loops in the / — v dia- 
gram; the extremity of a loop corresponding 
to a tangential direction in the galactic plane. 
The observed pattern is represented by the 
sum of the m = 2 system (two identical long 
arms with phase difference of 180°) and of the 
m — A system (four identical short arms each 
separated by 90° in phase) shifted by some 
phase angle from the first system. The ad- 
justed parameters are the pitch angles, the 
angle for the phase shift, the inner and the 
outer radii of each of the two systems. 

We used the catalog of H II regions of 
Kuchar & Clark (1997). The I - v dia- 
gram with the fitted loops, is shown in Fig.5. 
The rotation curve, that we adopted for the 
position-velocity transformation was derived 
from the interstellar gas data of Clemens 
(1985), reinterpreted in terms of i?o = 7.5 kpc. 
In this transformation we also took into ac- 
count a non-circular motion, due to the fact 
that when they form, H II regions have a sys- 
tematic velocity of about 10 km towards 
the center of the Galaxy, that keep them in or- 
bits which are close to the spiral arms (Bash, 
1981). Note that the choice of this parame- 
ter considerably helps in obtaining good fits 
of the observed locus of the inner HII regions. 
The loops are labeled a to z in the / — v dia- 
gram; the corresponding positions of the tan- 
gential points are indicated on the pattern in 
the galactic plane, in Fig. 6. 



We next discuss some of the features of 
the proposed structure. In the region /= 340° 
to 270°, which is the only region where clear 
arms were observed by GG, the structure re- 
sulting from our study closely resembles that 
of these authors, presenting about the same 
tangential directions. Remark that the longi- 
tudes of tangential directions arc not affected 
by a change of distance scale (GG used Rq = 
10 kpc). Therefore, in this range of longitude, 
it is almost impossible to distinguish between 
the 2+ 4 arms model that results from our 
study, and the empirical model of GG.This 
is specially true if we adopt a smooth func- 
tion to represent the potential, like we did in 
previous sections, since close potential min- 
ima tend to merge. A difference between GG 
and our work is that GG did not indicate the 
existence of a tangential direction at about 
Z=338°(our inner loop e), but obviously there 
are observed HII regions in that direction, 
at large negative velocities, that justify our 
model. In some other directions, the obser- 
vations favor our model as well. Remark for 
instance that there are concentrations of H II 
regions near labels a and / in Fig. 5. These 
are well explained by a spiral arm that pass 
very close to the Sun, seen almost at / ~ -|-90° 
and then at / ~ —90°, since the velocities 
are almost zero in these directions, for dis- 
tances that are not too large, according to 
the well known expression v oc sin{2l). The 
longitudes of the tangential directions indi- 
cate that it is an arm with small pitch angle 
(about 6°). On the contrary, the wide loop 
label b, can only be well reproduced with a 
larger pitch angle (12°). This emphazises the 
ncdd for arms with different pitch angles. The 
best fit of the / — v diagram of the observed 
HII regions that we obtained with a simple 2 
-I- 4 armed model, although not perfect, re- 
produces the main features of the diagram. 
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and in particular, the main tangential direc- 
tions. This was obtained with pitch angles 
6.6° for the 2-armcd component and 12° for 
the 4- armed component. 

Let us now discuss the theoretical / — v 
diagrams, that were computed by means of 
our particle simulations in Sec. 4 and are 
shown in Figs. 7 to 9. Figs. 7 and 8 are 
very similar to the observed diagram for H 
II regions, in many aspects. In particular, 
the observed loops like b, c, d, e which are 
interpreted in Fig. 5 as empirical locii of spi- 
ral arms, with corresponding tangential direc- 
tions to the arms, can be clearly seen in Fig. 
7 and 8 as well. A difference that appears 
between the theoretical I — v diagrams and 
the observed ones for H II regions is a strong 
concentration of points along the line from 
/ -60°, V ^ +100 km s-^ to / ^ +60°, v ^ 
— 100 km s^^. This is an expected result since 
the points corresponding to H I (recall that 
the particles simulate the H I clouds, see Sec 
4) are situated at large galactic radii, where 
H I is known to exist, whereas there is lack of 
H II regions at large radii. 

It is not surprising that the / — v dia- 
grams in Figs. 7 and 8 are so similar, since 
both models contain similar 2-arms compo- 
nent, and this component is prominent over 
a wider galactocentric region than the 4-arms 
component of the 2-|-4 arms model. However, 
we can point out a number of differences. For 
instance, the observed loop b, which is due to 
a 4-arms component, is correctly reproduced 
by the particles of the self-sustained model 
(Fig. 7) and does not appear in Fig.8. On 
the other hand, the arm that passes very close 
to the Sun, as discussed above (loops a and 
/) appear more clearly in the 2-arms model 
(Fig.8). We emphazise that the loops pre- 
sented in Figs. 7-9 are eye-guides for com- 



parison of theoretical diagrams with observed 
HII regions, but they are not perfect fits to the 
HII regions. In the HII regions diagram (Fig. 
5) we can see many objects between loops c 
and d, that are not well fitted by the loops. 
We can see many objects as well in this re- 
gion in Figure 8, the theoretical self-sustained 
model. In other words, the theoretical model 
is closer to reality than the "empirical" fit 
represented by the lines. Around longitude 
240° (or -120°), the theoretical I - v diagram 
of the self-sustained shows many objects with 
velocities of the order of 30 kms^^. These ob- 
jects seem to delineate a spiral arm that is 
not seen in the HII regions diagram. How- 
ever, a arm indeed exists at this position, as 
can be seen in the longitude- velocity diagram 
of IRAS sources given by Wouterlout et al. 
(1990). This arm is deficient in HII regions, 
probably due to the proximity of corotation. 

On the contrary, the 4 arms model (Fig. 9) 
show stronger differences with observations. 
The particles do not show loops a and b, but 
show a clear loop between loops c and d, that 
is not present in the HII regions diagram. 

In summary, the pure 2-arms model and 
the self-sustained 2-1-4 arms model produce 
theoretical / — v diagrams that are similar be- 
tween them and similar to that of observed 
HII regions. The only striking difference be- 
tween the theoretical and observed diagrams 
is an expected one, due to the fact that we 
are comparing objects that behave differently 
(HI gas and HII regions). If we look into the 
details, the 2+ 4 arms self-sustained model is 
favored. The choice in favor of the 2+4 arms 
model, compared to the pure 2-arms, is more 
strongly dictated by the study of Cepheid 
kinematics. 
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7. Conclusion 

In the present research a new approach to 
the problem of the galacic spiral structure was 
proposed in order to construct a more realis- 
tic picture like often seen in external galax- 
ies: co-existing of different spiral systems in a 
galaxy, arm bifurcation and their overlaping. 
Our theoretical considerations shows that su- 
perposition of self-sustained spiral wave har- 
monics could explain some of the above fea- 
tures, since different azimuthal wave harmon- 
ics have different pitch angles. 

In the framework of the simplest model of 
superposition of 2-|-4-armed spiral wave har- 
monics, wc analysed the best up-to-date data 
on stellar kinematics, which is the sample of 
Ccphcid stars, with proper motions and par- 
allaxes determined by HIPPARCOS. We ex- 
amined two models, the self-sustained and the 
bar-dominated waves. This study comple- 
ments the previous studies of pure 2-arms and 
pure 4-arms presented by MZ ad MZDMR. 
Of the four models, clearly the one that gives 
the best fit to the Cepheid kinematics is the 
self-sustained model, which is a superposition 
two arms with pitch angle about 6° and four 
arms with pitch angle about 12°. We per- 
formed N-particle simulations to make visible 
the structure of the potential derived from the 
Cepheid kinetics for the four models, and to 
construct the corresponding / — v diagrams. 

As an independent test of a 2+4 arms 
model, we performed a new analysis of the l—v 
diagram of the galactic Hll regions, which are 
the best tracers of the large-scale spiral struc- 
ture. We fitted the observed / — v diagram 
empirically with the locii of spiral arms, us- 
ing a 2+4 arms model. Coincidcntally, the 
best empirical fit was again found using pitch 
angles about 6° and 12°. We also compared 



the theoretical I — v diagrams from the parti- 
cle simulations with the locii of arms derived 
from Hll regions. Although the differences be- 
tween / — V diagrams of the theoretical mod- 
els (pure 2-arms, pure 4-arms, and 2 different 
models of 2+4 arms) are not striking, the 2- 
arms model and the self-sustained 2+4 arms 
models arc the ones that produce theoretical 
/ — V diagrams most similar to that of the Hll 
regions. 

Of all the arguments that we examined, the 
significantly better fit of the kinematics (the 
5"^ analysis) of Cepheids with the 2+4 armed, 
self-sustained model, is the most convincing 
one, but clearly the analysis of the Hll regions 
sample gives support to our interpretation. 
Although the Galaxy probably shows some 
deviations from any simplified model, the 2+4 
armed model with different pitch angles is the 
one that constitutes the best approach, being 
consistent with observations and with spiral 
waves theory. This is the simplest model that 
can be proposed, apart from pure harmonic 
modes. Our model does not exclude the exis- 
tence of higher modes, that are allowed to ex- 
ist near corotation, but these modes are prob- 
ably less significant than the first harmonics. 
It is interesting to remark that our model is 
able to reconcile the first model of Lin & Shu 
(1964), which has 2 arms with pitch angle 6° 
similar to our 2-arms harmonic, with the need 
to satisfy Kennicut's (1982) correlation be- 
tween pitch angle and maximum rotation ve- 
locity, which predicts a pitch angle of about 
14° for our Galaxy, not very different from 
that of our 4-arms harmonic. 

As a by-product of the study, our particle 
simulation shows that a deficiency of intertel- 
lar gas must occur near corotation. This ex- 
plains the gap in the HI ditribution observed 
by Kerr (1969) and by Burton (1976). This ef- 
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feet is also possibly related to the sharp min- 
imum in the rotation curve near the solar ra- 
dius discussed by Amaral et al. (1996), also 
seen in the curve by Honma & Kan-ya (1998). 
In external galaxies the ring-like gap in H I 
distribution may directly show the localiza- 
tion of the corotation circle. 

Our investigation of Cepheid kinematics 
points out that the Sun lies very close to 
the corotation circle. This conclusion is sup- 
ported by the solar closeness to the ring-hke 
region of H I deficiency in the Galaxy derived 
by Kerr (1969) and by Burton (1976), as well 
as by the position of Lindblad resonances (as 
discussed in Section 1) and by direct measure- 
ment of the pattern speed using open clusters 
(AL). 
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Figure captions 

Fig.l.- Surfaces of the A as a function of Xe2 
and for 3 values of pitch angle i2- 

Fig. 2.- The locus of minips2 and mimpsi on 
the galactic plane. The scale is indicated in 
kpc. As discussed in the text, this does not 
correspond necessarily to the visible struc- 
ture. 

Fig. 3.- The visible structure of the Galaxy de- 
rived for the best model (superposition of 2-1-4 
self-sustained wave harmonics) by means of 
particle-cloud simulation. 

Fig. 4.- The radial distribution of cloud con- 
centration < n > (averaged over a circle) for 
t = (dashed line) and t = 3.0 (solid line). 
The gap in the gas distribution is clearly seen 
near the corotation radius Re- The Sun is 
situated at Rq = 7.5 kpc. 

Fig. 5.- The observed / — v diagram for H II 
regions (from Kuchar & Clark, 1997) and the 
loops that we fitted empirically with a 2-1-4 
arms structure. 

Fig. 6.- The spiral structure of the Galaxy de- 
rived from H 11 data. 

Fig. 7.- The theoretical / — v diagram com- 
puted by means of particle-cloud simulation 
for the best model: superposition of 2-1-4 self- 
sustained wave harmonics.. The lines repre- 
sent the fit to observed HII regions, from Fig. 
5, for camparison. 

Fig. 8.- Same as in Fig. 7, but for the model of 
pure m — 2 wave harmonic (the parameters 
were taken from MZDMR).The lines repre- 
sent the fit to observed HII regions, from Fig. 
5, for camparison, and not a fit to the m = 2 
model. 

Fig. 9.- Same as in Fig. 8, but for the model of 
pure m = 4 wave harmonic (the parameters 
were taken from MZ)..The lines represent the 
fit to observed HII regions, from Fig. 5, for 
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camparison, and not a fit to the m — Amodel. 
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Table 1: Model parameters and their errors derived by means of statistical analysis. 



appro- 
ach 


/ km \ 
^ s kpc ' 


A 

f k7n \ 
^ s kpc ^ 


Rq^q 
/ kva \ 
^ s kpc^ ' 


(^) 


W0 

(^) 


n 


X02 

n 


/i?2 
(^) 


/i?2 
(^) 


K4I 

n 


X04 

n 


fm 
(^) 


hi 
(^) 


min 


self- 


26.3 


17.5 


9.9 


-8.8 


11.9 


6.1 


311. 


0.4 


-14.0 


12.0 


122. 


0.8 


-10.9 


187. 


sust. 


±1.3 


±0.8 


±1.9 


±1.0 


±1.1 


±0.4 


±11. 


±3.0 


±3.0 


±0.8 


±15. 


±3.3 


±2.9 




bar— 


25.7 


9.8 


-6.8 


-13.6 


9.3 


12.6 


184. 


12.5 


-19.4 


12.6 


8. 


6.6 


-10.1 


220. 


domin. 


±1.2 


±2.3 


±5.0 


±2.4 


±1.2 


±0.5 


±4. 


±3.8 


±4.4. 


±1.0 


±7. 


±2.4 


±2.1 
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